IMPERIAL /TP /06/SZ/Ol 



Evidence for an entropy bound from fundamentally 

discrete gravity 

D. Rideout, S. Zohren 

Institute for Theoretical Physics 
Imperial College 
London, SWT 2AZ, UK. 

E-Mail: |d.rideout@imperial.ac.uk[ |stefan.zohren@imperial.ac.uk| 



Abstract 

The various entropy bounds that exist in the literature suggest that spacetime is 
fundamentally discrete, and hint at an underlying relationship between geometry and 
"information" . The foundation of this relationship is yet to be uncovered, but should 
manifest itself in a theory of quantum gravity. We present a measure for the maximal 
entropy of spherically symmetric spacelike regions within the causal set approach to 
quantum gravity. In terms of the proposal, a bound for the entropy contained in this 
region can be derived from a counting of potential "degrees of freedom" associated 
to the Cauchy horizon of its future domain of dependence. For different spherically 
symmetric spacelike regions in Minkowski spacetime of arbitrary dimension, we show 
that this proposal leads, in the continuum approximation, to Susskind's well-known 
spherical entropy bound. 



1 Introduction 



In the history of general relativity there has been a long discussion regarding the thermo- 
dynamics of gravitational systems. One of the most famous examples is the Bekenstein- 
Hawking formula for black hole entropy El El E] , stating that the entropy of a black hole 
is given by a quarter of the area of the event horizon in Planckian units 

Sbh = J. (1) 

This bound is truly universal, meaning that it is independent of the characteristics of the 
matter system and can be derived without any knowledge of the actual microstates of the 
quantum statistical system. 

Along this line there have been several generalizations of this entropy bound. One is 
Susskind's spherical entropy hound [5^, stating that the upper bound for the entropy of the 
matter content of an arbitrary spherically symmetric spacelike region (of finite volume) is, 

>S'matter ^ (2) 

where A is the area of the boundary of the region.^ Even though this spacelike entropy 
bound cannot be generalized to arbitrary non-spherically symmetric spacelike regions, there 
exists a generalization in terms of light-sheets, namely the Bousso or covariant entropy 
bound jHl El- More precisely, let A{B) be the area of any {d — 2)-dimensional surface B, 
then the [d — 1) dimensional hypersurface L is called the light-sheet of i3 if L is generated 
by light rays which begin at B, extend orthogonally away from B and have everywhere 
non-negative expansion. The entropy flux through the light-sheet is then bounded by 

S{L) ^ (3) 



This entropy bound is also widely regarded as evidence for the holographic principle [H 
El Ej, stating that the maximum number of degrees of freedom carried by L is given by 
AiB)/A. 

These entropy bounds suggest that an underlying theory of quantum gravity should 
predict the bounds from a counting of microstates. This verification of the thermodynamic 
laws is an important consistency check for any approach to quantum gravity. Further, the 
finiteness of the entropy might already give some indications about the actual microstruc- 
ture of spacetime. There is a semi-classical argument that the description of a quantum 



^ There was an earlier proposed bound by Bekenstein [S] stating that the entropy of any weakly grav- 
itating matter system obeys S'mattor ^ 2-kER, where E is the energy of the matter system and R the 
circumferential radius of the smallest sphere that contains it. If one further assumes that this bound is 
valid for strongly gravitating matter systems, then gravitational stability in four dimensions implies that 
2E ^ R and hence S'mattor ^ 2itER < A/i. One sees that in four dimensions the Bekenstein bound is 
stronger then the Susskind bound, however, in d > 4 gravitational stability and the Bekenstein bound only 
imply that S'mattor ^ {d^2)A/8 7 . Hence, the geometrical Susskind bound is arguably more fundamental. 
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theory of gravity by a local quantum field theory in the continuum, in the absence of a 
high frequency cutoff, leads to infinitely many degrees of freedom in a finite region, and 
therefore to a divergence in the entropy of this region ■ The entropy bounds there- 

fore suggest that spacetime possesses a fundamental discreteness at scales of order of the 
Planck scale. Continuum physics would then have to emerge from this fundamental theory 
when making a continuum approximation at large scales. This suggests that to obtain a 
theory of quantum gravity one does not have to quantize the metric fields of the continuum 
geometries, but should rather find a quantum theory of the discrete structure underlying 
those continuum geometries [T!?j . 

In the following we introduce causal sets [HI [T31 Uni HZl HHI as an approach to funda- 
mentally discrete quantum gravity. We discuss the importance of causal structure within 
this approach, and give the correspondence between discrete causal sets and the emergent 
continuum manifolds in terms of faithful embeddings. We show how in this framework one 
can derive a notion of maximum entropy from a counting of potential horizon "degrees of 
freedom" of the fundamental theory reminiscent of former ideas in the context of black 
hole entropy f^lE]- Using this we formulate an entropy bound for spherically symmetric 
spacelike regions within the causal set approach. We then show that in the continuum 
limit, for different spherically symmetric spacelike regions in Minkowski spacetime of arbi- 
trary dimension, this leads to Susskind's spherical entropy bound. 



Conventions: 



Throughout this paper we work in Planck units with 

c=kB = GN = h=l, (4) 

where c, ks, Gat, h denote the speed of light, Boltzmann's constant, Newton's constant and 
Planck's constant respectively.^ In these units the Planck length is given by 

d-2 

I, = ' = 1, for O 3. (5) 



2 Causal sets: Fundamentally discrete gravity 

2.1 Discreteness and causal structure as first principles 

As mentioned in the introduction, causal set theory is an approach to fundamentally dis- 
crete quantum gravity. Besides taking fundamental discreteness as a first principle, the 
primacy of causal structure in the continuum is the main ingredient for causal sets. 

Causal order suggests itself as a fundamental principle for quantum gravity because of 
the enormous amount of topological and geometrical information which it contains [20] . It 

^It is perhaps more natural to set SttGn = 1, in which case ^ becomes Sbh = 2ttA, as is done in |19|. 
however in this work we keep with the more common convention Gat = 1. 
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has been shown that, given merely the causal relations of events in a spacetime manifold, 
one needs only the volume measure to recover the full geometry in the continuum fn\ • 
Furthermore, causality turned out to be a crucial ingredient for another approach to quan- 
tum gravity, namely Causal Dynamical Triangulations [23], which has led to considerable 
success e.g. over the Euclidean approach. 

In causal set theory, causal structure in the continuum is a reflection of an order rela- 
tion among causal set "elements". More precisely, a causal set is defined to be a locally 
finite partially ordered set C = (C, -<), namely a set C together with a relation -<, called 
"precedes" , which satisfy the following axioms: 

Transitivity: If x ^ y and y -< z then x -< z, Wx, y, z E C; 

Irrefiexivity: x -/^ x; 

Local Finiteness: For any pair of fixed elements x and z of C, the set of elements lying 
between x and z is finite, card{j/|a; -< y -< z} < oo, where cardX means the cardinality 
of the set X. 

Of these axioms, the first two say that C = {C, -<) is a partially ordered set or poset. 
The last expression, local finiteness, is important to recover the volume information of 
spacetime and is hence crucial for finding any continuum approximation of the causal set. 

In the following we will summarize some basic definitions in causal set theory of which 
most are "discrete" versions of the analogous concepts used to describe the causal structure 
of continuum spacetimes (cf. App.Ej). 

The past of an element x G C is the subset past(x) = {y E C \ y ^ x}. This corresponds 
to J~{x) in the continuum approximation. The past of a subset of C is the union of the 
pasts of its elements. The future of an element a; e C is the subset future(x) = G 
C\x -< y} which respectively corresponds to J^{x) in the continuum approximation. 

An important concept for the formulation of the entropy bound from causal set theory 
is that of a maximal element in a causal set C = (C, -<). This is an element which has no 
successors, i.e. an element x for which $y E C such that x y. The set of all maximal 
elements in C is denoted by max(C) = {x E C \iiy E C s.t. x -< y}. In analogy, a minimal 
element is one which has no ancestors, i.e. an element x for which $y E C such that y ^ x 
and the set of all minimal elements in C is denoted by min(C) = {x E C \$y E C s.t. y -< x}. 

2.2 Towards a continuum approximation 

Having given the precise definition of a causal set one might ask the questions: How can 
one actually formulate a theory of quantum gravity using causal sets, and how can causal 
sets lead to the known classical notion of smooth Lorentzian manifolds? 

The aim of causal set theory is to formulate a theory of quantum gravity using a 
sum-over-histories approach, where the single histories are causal sets. In order to recover 
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continuum physics in a semi-classical limit, some of the histories must be well approximated 
by Lorentzian manifolds. In the following we give a definition of what it means for a causal 
set C to be approximated by a Lorentzian spacetime {A4,g). 

Consider a strongly causal spacetime {A4,g). The map (p : C ^ A4 from a causal set 
C = (C, -<) into a spacetime {Ai,g) is called a conformal embedding ii x ~< y <^==^ G 
J^{(f){y)), Vx,y G C. Consider the Alexandrov neighborhood J^{p) fl J^{q), for every 
g G A^, which forms a basis for the manifold topology of {M.,g) if {M.,g) is strongly 
causal, which we assume throughout. The map (p is called a faithful embedding or sprinkling 
if it has the following property: The number of elements n mapped into an Alexandrov 
neighborhood is equal to its spacetime volume V, up to Poisson fluctuations. Thus, the 
probability of finding n elements in this region is given by the Poisson distribution 

P{n) = , (6) 

nl 

where p = Ij'^ is the density set by the fundamental length scale If in d dimensions. In 
other words, this means that in a spacetime region of volume V one finds on average 
= Vp elements of the casual set embedded into this region, and the fluctuations are 
typically of the order ^/N. Further, the properties of the Poisson distribution lead to local 
Lorentz invariance in the continuum spacetime [211 12S]- From the considerations above 
one assumes that the fundamental length scale will be of order of the Planck length. In 
Sec. m we will further predict an explicit value for the fundamental discreteness scale of 
spacetime which will be of the order of the Planck length. 

We say that a spacetime {A4,g) approximates a causal set C if there exists a faithful 
embedding of C into {Ai,g). This notion gives a correspondence between causal sets and 
continuum spacetimes. In terms of this correspondence continuum causal structure arises 
from the microscopic order relations of the causal set elements and the continuum volume 
measure of a region arises from counting the number of elements comprising this region. 
Nevertheless, one still has to prove the uniqueness of this correspondence (up to small 
fluctuations). In general the precise formulation of such a uniqueness proof is a difficult 
mathematical problem, and still remains a conjecture: the ^^Hauptvermutung" of the causal 
set approach. However, it has been proven to hold for the limiting case p i— > cx3 j2n] and 
certain progress has been made in the generalization to large but finite p |27j . 



3 An entropy bound from causal set theory 

In the previous section we have given a short introduction to the causal set approach to 
quantum gravity. Thereby, we gave a precise definition of the single histories of the sum- 
over-histories formulation at the microscopic scale. From a more phenomenological point 
of view one can raise the question of what conclusions one can draw from this for physics 
in the continuum. We have already seen, for example, how spacetime volume arises from 
a counting of fundamental spacetime elements. In the following we want to show how 
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Figure 1: Illustration of a spacelike liypersurface S, its boundary its future domain 

of dependence -D"'"(E) and its future Cauchy horizon H~^{T,). 



entropy bounds could arise from a counting of potential horizon "degrees of freedom" at 
the fundamental level, giving a microscopic origin for Susskind's spherical entropy bound. 

As already mentioned in the introduction, the spherical entropy bound states that 
the entropy of the matter content of a spherically symmetric spacelike region E (of finite 
volume) is bounded by a quarter of the area of the boundary of E in Planck units 

S < J, (7) 

in full units S ^ AkBC^/{4:GNh), where A = Vol(S(E)) is the area of the boundary of this 
region. 

Surprisingly, the maximum entropy in ((7j) can be determined without any knowledge 
of the microscopic properties of the thermodynamic system. A theory of quantum gravity 
however should be able to deduce ((Zj) purely from a counting of the fundamental degrees 
of freedom at the microscopic level. (One may wonder in what manner a counting of 
degrees of freedom measures the entropy of a system. In a discrete context "degrees of 
freedom" are generally finite, and a state counting can be expected to yield something 
proportional to the exponential of the number of degrees of freedom. Thus measuring the 
entropy as the logarithm of the number of states can be seen to be equivalent to counting 
the number of degrees of freedom of the system.) In the following we want to give a notion 
of those fundamental degrees of freedom in the context of causal set theory leading to the 
formulation of an entropy bound within this approach. 

Consider a spherically symmetric spacelike region of spacetime, potentially containing 
some matter system. We assume that the "back reaction" of the matter content upon the 
spacetime geometry can be neglected, so that different states of the matter system lead 
to the same spherically symmetric spacetime geometry. The entropy of that system must 
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eventually "flow out" of the region by passing over the boundary of its future domain of 
dependence. But because spacetime is fundamentally discrete, the amount of such entropy 
flux is bounded above by the number of discrete elements comprising this boundary. This 
is a fundamental limit imposed by discreteness on the amount of information flux which 
can emerge from the region, and therefore on the amount of entropy which it can contain.'^ 
Thus we have argued that an entropy bound arises from causal set quantum gravity. 

Proposal Consider a spherically symmetric spacelike hypersurface S of finite volume in a 
strongly causal spacetime M. of dimension d ^ 3. Denote the future domain of dependence 
of this hypersurface by D+(S) (c.f. App. El and Fig. Let C = {C, -<) he a causal set 
which can he faithfully emhedded into Then the maximum entropy contained in S 

is given hy the numher of maximal elements ofC, 

'S'max = card{max(C)}. (8) 

Claim This proposal leads to Susskind's entropy hound in the continuum approximation, 

Smax = (9) 

where A = Vol(i3(S)) is the area of the houndary of this region S, if the fundamental 
discreteness scale is fixed at a dimension dependent value to he calculated. 



4 Evidence for the claim 

In this section we provide analytical and numerical evidence for the claim that the entropy 
bound (jHl) leads to Susskind's bound in the continuum. In all discussed examples we con- 
sider {d— l)-dimensional spherically symmetric spacelike hypersurfaces S in d-dimensional 
Minkowski space M'^, and calculate the number of maximal elements in its domain of de- 
pendence D^CE). Since D+(E) is supposed to arise as a continuum approximation of a 
causal set C, we know from causal set kinematics that the elements of C are faithfully em- 
bedded into according to the Poisson distribution (jH)). Hence the expected number 
of maximal elements in D^(T.) is given by 

{n)=pf dx'^expl-pVol fj+(a;)nD+(S))| , (10) 

where again p is the fundamental density of spacetime. On the right hand side of (jlOp one 
integrates over all points x G Z?"'"(S), where every point is first weighted by the probability 
of finding an element of C embedded at this point and further weighted by the probability 
of not finding any other element of C embedded in J^{x) f]D"'"(S). Note that because of 

■^Those familiar with causal sets might worry that we are neglecting non-local effects when making this 
claim. However at some scale above the discreteness scale we must recover an effectively local dynamics, 
thus justifying our claim (perhaps up to some locality scale correction factor). 
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the fundamentally random nature of the discrete-continuum correspondence in causal sets 
we calculate the expected number of maximal elements in a sprinkling, even though the 
proposal (jHI) is phrased in terms of a fixed causal set. 



4.1 The ball in Minkowski spacetime 

In this section we want to calculate the expected number of maximal elements (n) in 
D"'"(S), where S is chosen to be a (ci — l)-dimensional ball Sd-i{R) with radius R in 
Minkowski spacetime M"^. Due to the spherical symmetry of the problem it is useful to 
introduce spherical coordinates x = (t,r,9i, ...,9d-2), where we choose the origin to be 
the futuremost event of Z^^(S). The volume element Vol(J^(a;) f]D^(S)) is equal to the 
volume of the Alexandrov neighborhood of x and 0, Vold(r) = Vol( J"*"(a;) f] J~(0)), where 
we denote proper time by r = \/t'^ — r^. Using the result for Volrf(r) as calculated in App. 
lEl and integrating out the spherical symmetry in (jiup one obtains a general expression for 
the expected number of maximal elements in D~^{Sd-i{R)), 



rR rt 



(n)=P ^^"2r / I drr'-^e-"^^^^'-^'^', (11) 
i [^-) Jo Jo 

where the dimension dependent constant is defined in App. ^ 

In the following we evaluate for various dimensions by analytical and numerical 
methods. 



4.1.1 2+1 dimensions 

In d = 3 dimensions one can explicitly evaluate (jlip . It is useful to express the re- 
sult in terms of the expected total number of spacetime elements faithfully embedded 
into D^{S2{R)), N = pV, where V = |-R^ is the volume of the domain of dependence 
{S2{R)) ■ The expected number of maximal elements is 

(n) = 8 (e-T - l) - 2NE^ (^^^ + AT (12) 

where En{x) is the exponential integral defined by 

/oo 
r'^e-'^'dt. (13) 

We also measured the expected number of maximal elements numerically by "sprin- 
kling" a causal set into a 2 -(-1- dimensional 2x2x1 square box, which contains D"'"(S'2(1)). 
By sprinkling we mean simply selecting Nc elements at random with uniform distribution 
within the box, and computing the causal relation between each pair from the Minkowski 
metric. For each of 100 trials i = 1 . . . 100, we deduce the set of elements which fall within 
D+(5'2(l)), compute its cardinality A^j, and count the number of such elements rii which are 
maximal within that region. From these we compute the sample mean and its error, and 
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Ball in 2+1 dimensions 
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Asymptotics 
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Figure 2: Shown is the expected number of maximal elements (n) as a function of the 
total number of elements in the domain of dependence of the unit disk in 2+1 dimen- 
sional Minkowski spacetime, on a logarithmic scale. Besides the analytical result and its 
asymptotics, data points with error bars from Monte-Carlo simulations are also shown. 



repeat this computation for a range of values of Nc- These computations were greatly fa- 
cilitated by utilizing causal set and Monte-Carlo toolkits within the Cactus computational 
framework [2H1- 

The plot of the expected number of maximal elements (n) as a function of for the 
unit disk is shown in Fig. |21 on a logarithmic scale. The agreement of analytical and 
numerical results justifies the the numerical methods. 

For a large number of elements ^ oo we can use the asymptotic expansion of the 
exponential integrals 

EJx) Qc 1 h... for Ixl ^ oo (14) 

X \ X / 

yielding 

(n) =4r j v^+... for oo, (15) 
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Figure 3: A snapshot from a simulation showing = 7806 spacetime elements forming the 
domain of dependence of the two-dimensional ball in 2+1 dimensional Minkowski spacetime 
(illustrated in red) and n = 126 maximal elements therein (illustrated in green). 



where + . . . are lower order terms in A^. The asymptotics are also displayed in Fig. El 
together with the full expression for the expected number of maximal elements. In terms 
of the density p and the radius R of 5*2 (-R) this result reads 



up to lower order corrections. It is important to see that (ji) ocy4/4, where A = 2'kR is the 
length of the boundary of 5*2 (-R). This is highly non-trivial, as one can see by looking at a 
snapshot of a numerical simulation (Fig. [21). There one observes that the maximal elements 
do not align along the one- dimensional boundary B{S2{R))- Instead they are distributed 
along a hyperbola close to the two-dimensional Cauchy horizon {S2{R)) , with a density 
of maximal elements which decreases with distance from the center. Hence the fact that 
the expected number of maximal elements is indeed proportional to the length A of the 
one-dimensional boundary B{S2{R)) for large A already gives very non-trivial evidence for 
the proposed entropy bound. 

To have the precise confirmation of the 0{1) constant in (n) ocA/A one has to choose 
the fundamental length scale // = p-^/^ to be // = 16v^/(35/6r(l/3)) ^ 4.41. This 
gives support to the belief that the fundamental discreteness scale of spacetime is of the 
order of the Planck length. Further, it gives a first determination for the coefficient in 2+1 
dimensions. 

Hence, using this value for the fundamental discreteness scale we have confirmed our 




(16) 



35/6r (1) ^' 4 
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claim, namely in the continuum approximation we have 



Smiix = {n) = (17) 

or S'max = kB^T^R/ {A^J KG / c^) in full units, when we set If = 16A^/(35/^r(l/3)). One 
now sees that the assumption of large but finite N used to obtain the asymptotic behavior 
is equivalent to saying that A is much larger than the Planck length. However, one can 
see that even for relatively small values of the length A, such as 10^ in Planck units, 
the approximation of (n) by its asymptotic expansion is already very accurate. For even 
smaller values of the length scale the Planck corrections become significant. However, as 
one can see in Fig. |2l the corrections always decrease the expected number of maximal 
elements, such that (n) never exceeds the bound (fTTj). 



4.1.2 3+1 dimensions 

Clearly from a physical point of view the evaluation of (jlip in 3+1 dimensions is the most 
interesting case. For c? = 4 one can write (jlip as follows, 

(n) = 27rp dt dzVW^e-P^"'' (18) 



Jo JO 

One can perform the first integration in ()18|) by using the following integral relation 

u 

x''-\u - x)^-^e^^"dx = 5(/i, z/)m'^+^-^ X 



X A ( .... ^i±il±l. .... ^ + + " - ^ ^pA . (19) 

\n n n n n n J 

for 3?(/i) > 0, 3?(z/) > and n = 2, 3, where B{fi, v) denotes Euler's beta function and 
pFq{ai, ...,ap;bi, ...,bq; z) is the generalized hypergeometric function defined through 

OO ^ T-rp / \ 

Ma„ ap- h, bq- z) = J2 iTTrHrr' ^20) 

fc=o \.\.j=iW)k 

and {a)n = r(a + n)/r(a) are the usual Pochhammer polynomials. The second integration, 
namely the one of the generalized hypergeometric function, can be obtained by use of the 
following relation [22] 

/z°' 
z"~^ pFg{ai, ...,ap; bi, ...,bg; z)dz = — p+iFg+i(ai, ...,ap,a;bi, ...,bg,a + 1; z). (21) 

Expressed in terms of the number of causal set elements sprinkled into D"*" (6*3(7?)), = pV, 
where the volume of D~^{Ss{R)) is given hj V = f-R"^, the final result reads 

W = iV3i^3Q,l,l4^,2;-iiv). (22) 
11 
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Figure 4: Shown is the expected number of maximal elements (n) as a function of the total 
number of elements in the domain of dependence of the 3-dimensional unit ball in 3+1 
dimensional Minkowski spacetime on a logarithmic scale. The plot shows the analytical 
result, its asymptotics, and numerical results from Monte-Carlo simulation. 



For a large number of elements we can use the asymptotic expansion of the gener- 
alized hypergeometric functions (cf. App. [H)), yielding the asymptotic expression for the 
expected number of maximal elements, 

{n) = 3V2^N^ + ... for A^ ^ oo. (23) 

In terms of the fundamental density of spacetime p and the radius R of 6*3 (i?) this result 
translates into 

(n) = ^pi^. (24) 

As in the previous case the number of maximal elements follows the right scaling according 
to the entropy bound, namely (n) oc A/4, where the area of the boundary of 5*3 (i?) is given 
hj A = AttR^. The factor of proportionality is a 0{1) constant as in the previous case, 
supporting the assumption that the fundamental length scale // is proportional to the 
Planck length Ip. More precisely, to have an exact agreement with the Susskind bound 
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the fundamental length scale in four dimensions will be given by // = in Planck units. 
Using this value for the fundamental length scale the asymptotic expansion for the expected 
number of maximal elements reads 

5'max = {n) = (25) 

or S'max = 4:7TR'^kBC^/{4:hG) in full units. At first sight one might feel uncomfortable in 
absorbing the order one constant into the fundamental length scale to exactly confirm 
the Susskind bound. However, as discussed in the previous case, the scaling (n) ocA/A is 
nontrivial and already serves as a confirmation of the bound. Taking the phenomenological 
law Smax = A/4: (spherical entropy bound) as "data" gives the fundamental length scale 
in four dimensions to be 

// = ~ 1.57 (in four dimensions). (26) 

In comparison to the previous case of 2+1 dimensions one observes that the fundamental 
discreteness scale // depends on the dimension. In Sec. I4.1.il we will derive an expression 
for the factor for arbitrary (even) dimensions, showing that this factor tends exactly to 
one a.s d —>■ oo. It is important to check that // is universal for all cases with the same 
spacetime dimension. In Sec. 14.21 we will provide evidence that this is indeed the case. 

At this point it is interesting to note that a similar method of fixing the fundamental 
discreteness scale to obtain the right factor of proportionality in the entropy bound is 
followed in loop quantum gravity in the context of black hole entropy (cf. ISO])- There 
one fixes the Immirzi parameter which can be regarded as a measure for the discreteness 
scale (through its relation to the lowest eigenvalue of the area operator) to obtain the right 
factor of a quarter in the black hole entropy. However, since there are several ambiguities 
in the relation between the Immirzi parameter and the fundamental discreteness scale in 
the sense one uses it in causal set theory, it is hard to compare the numerical values coming 
from the two derivations. 



4.1.3 4+1 dimensions 

In (i = 5 dimensions one can also evaluate the integral in (jllll in a similar way to the 
calculation in 2+1 dimensions. As in the previous cases the result is expressed in terms of 
the number of causal set elements sprinkled into D^{S4{R)), N = pV, with the volume of 
D+{S4{R)) given by ^ = ^R^. The final result reads 

+ 5(2Ar)tr (^0 - 30{2N)T>T (^0 ^ . (27) 

The plot of (n) as a function of N is shown in Fig. on a logarithmic scale as well as 
the numerical results obtained from Monte-Carlo simulation and the asymptotic behavior. 
As in the previous cases there is agreement between analytical and numerical results. 
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Figure 5: Shown is the plot of the number of maximal elements (n) as a function of the 
total number of elements in the domain of dependence of the 4-dimensional ball in 4+1 
dimensional Minkowski spacetime on a logarithmic scale, its asymptotics, and numerical 
results. 



For large one can expand (j27|) using the asymptotic expansion for the exponential 
integral (fT^ . yielding 

i^n) = —v(-\{2N)l + ... for N ^ oo. (28) 



Re-expressed in terms of the density p and the radius R of S^i^R) this result reads 



327ri/5 2 3 27r2i?3 

As for the lower dimensional cases this result gives the correct behavior [n) ocy4/4, where 
A = 2tt'^R^ is the volume of the boundary of S^lR), providing further evidence for the 
claim. Further, we can use the coefficient in (j29|) to fix the fundamental length scale in 
4+1 dimensions // = p'^/^ to be // = 2^^/^7r^/^y{5^/^^{5 + ^5)^/^1(8/5)^/3) ^ ^ 22 in 
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Planck units. Using this fundamental length scale the number of maximal elements ()29|) 
reads 

^Jmax = {n) = , (30) 

or in full units S'^ax = 27r'^ R'^kB/{4:(hGN/c^Y^'^), which confirms our claim. 



4.1.4 Generalizations to higher dimensions 



In the previous sections we have seen that causal set theory can provide a fundamental 
explanation for Susskind's entropy bound for the spacelike hypersurface Sd-i{R) in 2+1, 
3+1 and 4+1- dimensional Minkowski spacetime. In this section we want to generalize these 
results to arbitrary dimensions. For the case of odd spacetime dimension, it turns out that 
one can easily calculate the expected number of maximal elements, but one cannot write 
the results in a closed form. However, for even dimensions one can find a closed expression. 

As in the previous cases we will express the result of the expected number of maximal 
elements in terms of iV = pV, where V is the volume of the domain of dependence of 
Sd-i{R), given by 



d-l 

TT 2 



(31) 



Using the integration relations (fTIHl and ()21|1 one can integrate (fTTjl for even dimensions 
yielding 

'2 4 d ^ d + 1 d + 3 2d -I 
d' d 



In) = N d,,Fd 



+1 
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d 



2; -2^-'^A^ 



(32) 



Note that this result is also valid for the case of d = 2 dimensions, however it is not related 
to any entropy of the system, and is thus excluded from the proposal. 

For large one can use the asymptotic expansion of the generalized hypergeometric 
function (App. EI) to derive the asymptotics of the number of maximal elements for 
arbitrary even dimensions, yielding 

2—TT{d-l) 
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for A^ —>■ oo. 



From this result one can see that the number of maximal elements scales like (n) 
where A is the volume of the boundary of Sd-i{R), i.e. 
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'A/ A, 

(34) 



Further, we obtain 
scale. 



A = Yo\iBiSd-,iR))) 
(n) = A/4 precisely for the following value of the fundamental length 
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Figure 6: Shown is the plot of the fundamental length scale If in Planck units as a function 
of the dimensions d for even and odd dimensions and the analytic continuation. 



The analytic continuation of (j35|) as a function of the dimension is shown in Fig. Altogether 
with the explicit values for 2+1 and 4+1 dimensions as determined in the previous sections. 
One observes that the expression agrees with these values. This suggests that 
also holds for arbitrary odd dimensions. For d = 2 the value of diverges, since the 
Planck length lp = (IiGn/c^)^'^"'^^^'^ is not well defined in two dimensions. This also reflects 
the fact that the entropy bound only holds for dimensions d^3. For all values d^3 the 
fundamental length scale is of order of the Planck length. Interestingly, for large values 
d^oo expression (jH^ tends exactly to one. It is also less than one for values d>do, where 
do -10. 

4.2 Generalizations to different spatial hypersurfaces 

In the previous section we have derived the Susskind bound for the case where the spacelike 
hypersurface was chosen to be a (d — l)-dimensional ball in d dimensional Minkowski space- 
time. Further, from this we determined the fundamental discreteness scale of spacetime. 
However, it is important to prove that the fundamental discreteness scale so determined 
yields the same entropy bound for all spacelike hypersurfaces of a certain dimension. 
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Figure 7: Illustration of the different hyperbolic spherically symmetric spatial hypersurfaces 
S parameterized by a, and the domain of dependence 



In this section we show that the claim also holds for spacelike hypersurfaces in Minkowski 
spacetime different from those discussed in the previous section. We investigate hyperbolic 
spherically symmetric spacelike hypersurfaces given by 

t = -^r"^ + a?{R^ -r"^), ^ a ^ 1, (36) 

as shown in Fig. [7| together with its domain of dependence. For the special case of a = 1 the 
spacelike hypersurface given by (jHHjl is equivalent to the {d — l)-dimensional ball Sd-i{R) 
for which we have determined an analytic expression for the expected number of maximal 
elements (n) in the previous section. For other values < a < 1 one cannot determine the 
number of maximal elements analytically. In the following we will investigate this problem 
numerically for the physically most important cases of 2+1 and 3+1 dimensions. 



4.2.1 2+1 dimensions 

For d = 3 dimensions we use Monte-Carlo methods to numerically obtain the number of 
maximal elements in the domain of dependence of the spacelike hypersurfaces S defined 
by p6|) for different values of a as a function total number of elements in the domain of 
dependence. Since all these spacelike hypersurfaces S have the same boundary it is 

useful to express (n) as a function of the length of the boundary A = 2ttR. One can do this 
by using that N = pV, where the volume of the domain of dependence of S is given by 

27r o , , 

V = -^^R'. 37 
3 1 + a 
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Figure 8: Shown is the plot of the expected number of maximal elements (n) in the domain 
of dependence as a function of the length of the boundary A for different hyperbolic spher- 
ical symmetric spacelike hypersurfaces parameterized by a in 2+1 dimensional Minkowski 
spacetime. One sees that all functions approach the same asymptotic {n) = A/ A. 



Further, we use the value for the fundamental density of spacetime as obtained in Sec. 
HXn i.e. p=(35/2r(l/3)3)/(81927r). 

The results of the simulations are summarized in Fig. |H1 Shown is the expected number 
of maximal elements (n) as a function of the boundary length A. For the special case 
of a = 1 in ()36p we know that the result is given by ()12j). where N is replaced by = 
A^/(2^^r(2/3)^). The asymptotics of this analytic result is given by (n) = A/ A as shown 
earlier. For different values of < a < 1 the simulations show that even though the number 
of maximal elements as a function of A differs for small values A < 10^ in Planck units, in 
the expansion for large A, corresponding to the continuum approximation, all functions {n) 
for different a enter the same asymptotic expansion (n) = A/ 4, yielding Susskind's entropy 
bound. In addition, it shows that the prediction of the fundamental discreteness scale is 
universal in 2+1 dimensions, at least for all investigated cases of spacelike hypersurfaces 
in Minkowski spacetime. A generalization to examples in curved spacetime has still to be 
shown and will be investigated in future work. 
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3+1 dimensions 
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10 100 1000 

area of the boundary A 

Figure 9: Shown is the plot of the expected number of maximal elements (n) as a function 
of the length of the boundary A for different hyperbolic spherical symmetric spacelike 
hypersurfaces in 3+1 dimensional Minkowski spacetime parameterized by a. One sees that 
all functions approach the same asymptotics (n) = A/A. 



4.2.2 3+1 dimensions 

As in the case of 2+1 dimensions we use Monte-Carlo methods to numerically obtain the 
number of maximal elements in the domain of dependence of the spacelike hypersurfaces 
S defined by fl^ for d = 4 and different values of a. Again, the result is expressed as a 
function of the area of the boundary of S, i.e. A = 4:TtR^. This can be done by using the 
relation = pV and noticing that the volume of the domain of dependence of E is given 
by 



V = R'' 



TT 



TT 



6(1 



a?) 2 



Vl -a2(2 - Sa^) + 3a^log 



1 + v^ 



(38) 



where we use p= 1/6 for the value of the fundamental density of four dimensional spacetime. 

The numerical and analytical results are shown in Fig. IHl Displayed is the expected 
number of maximal elements (n) as a function of the area A of the boundary. For the special 
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case of a = l the analytical result was given by ()22|). with replaced by A^ = A^/(27rl2^). 
The simulations show that for different values of < a < 1 all functions (n) enter the same 
asymptotics, (n) = A/A, in the continuum limit, giving further evidence for the claim. 
For very small values of a the expected number of maximal elements enters the asymptotic 
regime only very slowly. However, the upper value for A of the simulations is still very small 
(~ 10~^'^m). As in the lower dimensional case the simulations show that the prediction 
for the value of the fundamental discreteness scale of four-dimensional spacetime, namely 
// = v^, is a universal quantity for this dimension, at least for all investigated cases of 
spherically symmetric spacelike hypersurfaces in Minkowski spacetime. 

5 Conclusion 

In this work we motivated causal sets as a candidate for a theory of fundamentally dis- 
crete quantum gravity. Causal structure and discreteness are implemented as fundamental 
ingredients in this approach and local Lorentz invariance follows as a consequence. In 
contrast to attempts to quantize general relativity by quantizing the metric fields, such 
as done in loop quantum gravity, causal set theory aims to formulate quantum gravity as 
a sum-over-histories, where the single histories are not smooth Lorentzian manifolds but 
rather ordered sets of spacetime elements. Continuum physics is then thought to emerge 
from this discrete structure in a continuum approximation. 

Using the kinematics of this framework we argued for a bound of the entropy in a 
spherically symmetric spacelike region from a counting of potential horizon "degrees of 
freedom" of the fundamental theory, namely the maximal elements of the future domain 
of dependence of the region. It was then shown that, for different spherically symmetric 
spacelike regions in Minkowski spacetime of arbitrary dimension, this leads to Susskind's 
spherical entropy bound. The evidence was given in terms of analytical results for spatial 
[d — l)-dimensional balls in d-dimensional Minkowski spacetime, and in terms of numer- 
ical results obtained from Monte-Carlo simulations for the case of hyperbolic spherically 
symmetric hypersurfaces in Minkowski spacetime. Agreement with the factor of propor- 
tionality in Susskind's entropy bound required the fundamental discreteness scale of four- 
dimensional spacetime to be // = 1.57 in Planck units. This is a "prediction" for 
the value of the fundamental discreteness scale of four dimensional spacetime within the 
causal set approach. (The first attempt of a determination of the discreteness scale came 
from computing the entropy of dimensionally reduced black holes, both Schwarzschild and 
spherically symmetric gravitational collapse, as described in |1JJ.) Even though this value 
varied with dimension, it was shown that, for all investigated cases, its value did not de- 
pend on the choice of spherically symmetric region for a given dimension. However, the 
evidence provided so far was given only in terms of examples in flat Minkowski spacetime. 
Clearly this is a very restricted class and further evidence should be provided in differ- 
ent spacetimes. Spherically symmetric spacelike regions in curved spacetime such as in 
Friedmann-Robertson- Walker cosmology will be investigated in future work. Another step 
would be to formulate the proposal intrinsically in terms of order invariants, without mak- 
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ing explicit reference to the continuum. Such a formulation may be fruitful in providing a 
fundamental understanding of Bousso's covariant entropy bound. 
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A Basic definitions regarding the causal structure 

In this appendix we state some of the basic concepts regarding the causal structure of 
continuum spacetimes. For further references the reader is referred to |SI1E21; where we 
follow the conventions of 

Let {A4,g) be a time orientable spacetime. We define a differentiable curve X{t) to be 
a future directed timelike curve if at each point p G A the tangent t° is a future directed 
timelike vector. Further, X(t) is called a future directed causal curve if at each pE X the 
tangent is a future directed timelike or null vector. 

Using this definition one can define the chronological future (past) and causal future 
(past) of a spacetime event or a set of spacetime events. The chronological future of p^Ai, 
denoted by I~^{p), is defined as 

I~^{p) = {q E Ai \ 3 future directed timelike curve X(t) s.t. A(0) = p and A(l) = q} (39) 

For any subset T, G A4 we thus define 



In analogy to I^{p) and /"""(S) we can also define the chronological past I~{p) and /~(S) 
by simply replacing "future" by "past" in (jHTIjl . 

In analogy to the chronological future we can also define the causal future of a point 
pG A^, denoted by J^ij)), by replacing "future directed timelike curve" by "future directed 
causal curve" in ()39|) . The definition of J^(S), J^{p) and J^(S) then follow accordingly. 

Very important for the formulation of the entropy bound stated in Sec. 121 is the concept 
of the domain of dependence of a spacetime region. 

Let E be any spacelike hypersurface of M.. We define the future domain of dependence 
of S, denoted by D"'"(S), by 

D^(T?) = {p E M. |every past inextendible causal curve through p intersects S} . (41) 




(40) 



pes 
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Analogously, one can define the past domain of dependence of S, denoted by D^{T,), by 
simply replacing "past" by "future" in (jUJ. 

Using the definitions above one can then define the future and past Cauchy horizon. 
Let S be any spacelike hypersurface of Ad. The future Cauchy horizon of S, denoted by 
is defined as 

iJ+(S) = D^\/-(Z}+(S)), (42) 

where D+(E) is the closure of The past Cauchy horizon, H^{I1), can be defined in 

analogy to (@2I)- The different concepts defined in this subsection are illustrated in Fig. ^ 



B Volume of the causal region between two events 

In this appendix we calculate the volume of the intersection of the causal future of an 
event p with the causal past of another event q G J^{p) in Minkowski spacetime M'', i.e. 
the Alexandrov region J^{p) {^J^{q). This volume depends only on the proper time r 
between p and q and hence we set Vol(i(r) = Vol( J+(p) J^{q)). 

Note that the volume of a d-dimensional ball with radius r, S'd(r), is given by 

^oM^^^^:C,r', „.th C.^^. (43) 



Using this we obtain 

Vol,(r) = 2 r dtYo\{Sd-i{t)) 
Jo 



d-l 



■.Ddr\ with Drf=-£-i-. (44) 



— ; -jTi—T ='■ J^dT , witn Ud = , , 

2'i-^dT{^) 2<^-^d 



2 

Specifically we get Vol2(r) = \t^, Vol3(r) = f^r^ and Vol4(r) = fj\ 



C Asymptotics of the generalized hypergeometric func- 
tion 

In this appendix we obtain the asymptotic expansion of the generalized hypergeometric 
function as given in (jH^ to first order. 

We use the general expression for the asymptotic expansion of the generalized hyper- 
geometric function pFq with p = q which can be derived from jSH] 

pFp(ai, a,; 6i, h,- z) = [f{ ^ j e^z^ + Q (^-^^ + 

^ (n^)E n-J(.-..) '->-0-Ki))' 
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where 

p 

7 = ^ - 6fc and V j. (46) 

k=l 

In the special case of flH^ we cannot straightforwardly apply (jlSjl . since ap_i = ap. However, 
since we are only interested in the first order behavior and ai < aj with ai 7^ aj for all 
2^j^p, we can obtain the asymptotic expansion of (j32p to first order, yielding 

/2 4 d d + 1 d + 3 2d-l , 



2^vr 



sin(f)r(^) 



AT-d +c»(Ar-d), AT ^00. (47) 
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